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discuss the connection between our results and the standard description of decoher- 
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(even when the correlations occur) reduced dynamics can be written in a manageable 
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I. INTRODUCTION 



Recently, the connection between a problem of decoherence and the Riccati operator 
equation was established (for details see^ and Ref. therein). Moreover, it was shown that a 
wide class of a time-dependent quantum systems, precisely the ones that describe a qubit 
Q interacting with the environment E and defined by the following Hamiltonian 

HQE{t, P) = Hgit, f3)^lE + lQ^HE + /(aa) ® V, (1) 
where the Hamiltonian HQ{t,l3) of the qubit alone is give by 



Hqit, (3) = (3a3 + a (ai cos{ujt) + a2 sm{ujt)) , (2) 

where a, /3 G M, can be effectively simplified to the time-independent problem. Namely the 
one governed by the Hamiltonian Hqe{0, (3) = Hqe- The connection between the reduced 
dynamics ' of those models can be summarized in the following formula 



(3) 



where (3 := (3 — co/2 plays the role of effective parameter and Vt := diag(e~*^*/^, e*^*/^) is 
the unitary (similarity) transformation. Here, pt is the solution (reduced dynamics) of the 
system identified with the time-dependent Hamiltonian HQE{t, [3) and pt{[3) representing the 
reduced dynamics of the model described by the Hamiltonian Hqe- An explicit dependence 
pt of (3 was omitted. The meaning of the symbols we used in the equations (1) and (2) is 
usual. Furthermore, it was found that this time-independent problem can be solved in two 
ways. One could either resolve for X the Riccati operator equation: 



aX"^ + X{H^ + 13)- {H^ - P)X - a = 0, 



(4) 



where H± := He ± V, or solve the following Schrodinger equation (we work with the units 
h=l) 



t'^t = Ht^t, Ht 



He zlVp 
ZtVp He 



(5) 
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with initial condition = ^- In Eq. (5) = [^i, 4>tf e H ® H, zt = e"'^"* and Vfj = 
V + (3Ie- In the description above it was assumed that "H is a separable Hilbert space 
(possibly infinite-dimensional) and He and V are the Hermitian operators acting on it. One 
could recall that in the current paper we assume that the Riccati equation Rh[X] = 0, 
where 



Rh[X] = XBX + XA-CX - B\ 
is associated with the following Hamiltonian 



(6) 



H 



A B 
5t C 



(7) 



In turn if the solution of the equation Rh[X] = exists, it may be used to diagonalize 
operator matrix H. Following equality holds true 



SJ^HSx 



A + BX 
A-{XBy 



(8) 



where 5*^^ stands for the inversed operator matrix to 



S 



X 



(9) 



I^, -xt 

X Ie 

The Riccati Eq. (4) is associated with the Hamiltonian Heq- On the other hand the Riccati 
equation related to the Hamiltonian reads 



X{z;V^)X + XHe - HeX - ztVp = 0. (10) 

One should keep in mind that the given operator A acting on the space ® "H may be 
thought of as a block operator matrix, since the following isomorphism holds C^^'H = 'H(B'H. 

It is worth mentioning that in spite of the fact that connection between operator matrices 
Hf and Hqe is well defined and the solution of the Eq. (10) can be easily obtained (indeed, 
it is given by Xt = ZiIe ) the solution of the Eq. (4) is still missing. Please note that 
although Xt is known it cannot be effectively used to resolve Eq. (5) because of its explicit 
time dependence. Nevertheless it allows us to diagonalize the Hamiltonian Ht, namely 
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H+ + (3Ie 

i/_ - ^Ie 



where the unitary matrix Szt is defined as 



V2 




We want to emphasize that the Eq. (10) and its shght modification 



(11) 



(12) 



zlXVpX + XHe - HeX - ztVp = 0, 



(13) 



used by the author in previous manuscript on the subject^, are equivalent only if the solution 
is assumed to be a linear operator. This seems to be justified, especially if one expects X 
to represent an observable. However, it does not need to be true. Thus in this manuscript 
we will not restrict the analysis to the linear operators only. 

The purpose of this manuscript is two fold. Firstly, we show that if the environment 
is {i.e., the operators He and V are) invariant under antylinear transformation (Sec. II), 
then this antylinear operator is the solution of the Riccati equation (10). We also indicate 
that mentioned symmetry may be used to diagonalize the Hamiltonian related with Eq. (4), 
although it does not solve this equation. Next, in Sec. (Ill) we argue that the solution 
we obtained can not be applied to solve Eq. (5). The problem occurs because standard 
methods provided by the theory of the differential equations demand that the operator is 
not antylinear. The reduced dynamics of the system under consideration is given in Sec. (V). 

Secondly, in Section (IV) we study the case when [He, V] = and we give the exact 
solution of the Riccati equation (4) for this situation. This is direct generalization of the 
recently found solution for the particular operator defining the spin-environment. Finally, 
in Sec. (VI) we compare our method with the standard approach based on the operator 
sum representation. We also discuss the possibility of obtaining an operator (or Kraus) sum 
representation in the case when the correlations between the system and its environment 
are present initially. 



II. SYMMETRY AND THE SOLUTION TO THE RICCATI EQUATION 



Let Ti and T2 be an antilinear symmetry for He and V, such that Tj^ = I^, (or Ti = r~^) 
for i = 1,2, respectively. By definition, it means that [He, ti] = and [V, T2] = 0. Since the 
symmetry operator is an involution i.e., = I^;, thus for a given operator A the condition 
[t. A] = is equivalent to the equality tAt~^ = A. Operators that fulfill last equality are 
invariant under the action of r and are called r— symmetric. Furthermore, the statement 
that an operator r acting on "H is an antilinear has the following meaning 

T{aip + b(j)) = a*T'ip + b*T(f), (14) 

for every tp, (p and a, 6 G C. We wish to emphasize that in the finite-dimensional case the 
existence of the aforementioned symmetry is ensured by the theorem of Ali Mostafazadeh'^ 
which states that every diagonalazable pseudo-Hermitian (in particular the Hermitian) op- 
erator H with the discrete spectrum has an antilinear and anti-hermitian symmetry r. 
Moreover, this symmetry is an involution, i.e., = Ie, if H is Hermitian operator, i.e., 
H = HK The proof of this theorem provides the explicit construction of r. Let us addition- 
ally assume that ti = T2 := r, which means that the operators He and V posses a common 
symmetry r. 

If one allow X to be an antilinear operator, then equation (10), can be rewritten in the 
following form 

Rh, [X] = zt {XVpX - Vp) + [X, He]. (15) 

From this equation we can readily see that R[t\ = 0, i.e., the symmetry generator r of the 
E system is the solution we were seeking for. In particular, if the Hamiltonians He and V 
are (or roughly speaking the E system is) PT— symmetric (or T-symmetrix) where P and T 
stand for the parity and the time- reversal operators, respectively, then X = PT {X = T), 
i.e., the PT (T) operator solves Eq. (10). This is very interesting and rather unexpected 
result that the famous PT (T) symmetry is the solution of the Riccati equation we study. 

For instance it can be easily proven that if both He and V are symmetric operators, that 
is to say H^ = He and V'^ = V, where by "T" we denoted the operation of transposition, 
then it implies that RHt[K] = 0, where K is complex conjugate operator (see below). In 
order to do that let us define 
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Kij = ^p\ e (16) 

where by "*" stands for the standard complex conjugation operation. The operator above 
possesses the following properties: 

a) K = {i.e., it is Hermitian operator), 

b) KK^ = Ie (i.e., it is unitary operator), 

c) = I^; (i.e., it is an involution). 

Moreover, K is antilinear operator. The listed properties follow immediately from the defi- 
nition (16). Note that for any Hermitian operator (matrix) A the condition = A means 
that A is i^— symmetric, i.e., [A, K] = 0. To see this observe that K transforms any operator 
A in accord with the following rule: 



KAK^ = A*. (17) 

In addition, if one assumes that A is Hermitian, i.e., A = and uses properties a) then 
from the equation above we have 



KAK = A^, (18) 

where we used the fact that A^ = (A^)*. Clearly, for the symmetric matrix KAK = A or 
[A, K] = 0. Other way to see that RHt[K] = is to rewrite R[K] as 



RhAK] = ZtiV^ - V^) + (Hi - He)K. 



(19) 



Evidently, for the symmetric operators the right side of the Eq. (19) vanishes. 

At the end of this section we show how to diagonalize the operator Hqe using the sym- 
metry r. Firstly, note that if one introduces the unitary operator f/ in a way that 



V2 



Ie He 
iIe Ie 




(20) 
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then WHqeU = H and 



H 



R 



(21) 



where Vap := als — iVp- Since rVapT = Vj^, thus from Eq. (21) we see that Rh[t] = 0, i.e., 
the symmetry r is the solution of the Riccati equation associated with the Hamihonian (21). 
Therefore, the matrix Sr '.= USr diagonahzes the Hamiltonian Hqe. To be specific, the 
following equation holds true 



Keep in mind that the operator Va^ is not Hermitian, while the operator VapT is. Interest- 
ingly, we diagonalized the block operator matrix Hqe without directly resolving the Riccati 
algebraic equation associated with it. One may ask if this is possible in general. To be spe- 
cific if there exists a matrix T such that it transforms a given operator matrix A to another 
one A viz A = T^^AT and the Riccati equation -R^[X] = related to A has the solution 
that is easy to be found. If that would be the case, than one may construct similarity 
operator matrix, namely Sx^T that diagonahzes A. We saw that if A = Hqe, then T = U 
(Eq. (20)) and A = H (Eq. (21)), thus the procedure is working. Therefore, this strategy is 
worth to be studied. However, we will not focus on this subject herein. 

III. THE PROBLEM WITH THE SCHRODINGER EQUATION 

The results of the paper^ and the last section show that RH^[zt] = and i?i^Jr], where 
RHt[X] stands for the left side of the Eq. (10). Note that, the second solution {X = r) we 
obtained, in contrast to the first one {Xt = Zt) is time independent. Therefore, one may 
think that the second solution has the advantage because it allows one to construct the 
operator matrix Sr that diagonahzes Hamiltonian (5) and do not depend on time. Indeed 
this is the case however one serious drawback arises. To see this clearly, let us give an 
explicit form of the matrix Sr. In agreement with Eq. (9) it takes the form (compare this 
with Eq. (12)) 




Hr + V^pT 




(22) 



Ob 
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s -1. 



Ie -r 
r Ie 



(23) 



Since the symmetry r is an involution operator (r^ = I^), thus Sr is the unitary matrix 
{St-SI = Ie)- We see that the similarity transformation (23) is neither linear nor antilinear 
operator. In particular, one may easily verify that Sri = iS^. All this difficulties are 
direct consequence of the fact that the operator r is antilinear. As a result, we cannot use 
the standard procedure that allows us to solve linear, differential equation to resolve the 
Schrodinger equation z|\E't) = Hfl'^t). Indeed, because of the presence of the factor i on the 
left side of the Schrodinger equation we cannot apply the \^t) = Srl'^t) transformation to 
reduce it to the following form i\^t) = Hf\^t)i where Hf stands for the diagonal form of Ht- 
Another technique is needed to resolve this difficulties. Unfortunately, at the present time 
there is none. Of course, the same problems occur when one tries to resolve the equation 
i\i!t) = HQE\^t) using the Eq. (22). 

Notice that according to procedure explained in the paper^ the Ht given by Eq. (5) has 
the following block-diagonalization (compare with Eq. (11)) 



S\HtSr 



He + ztVpT 








= H^ 



t ■ 



(24) 



He - ZtVpT 

Comparing equations (11) and (24) one can learn that the solution Xt = Zt leads to the time- 
independent diagonal form H"^ of the operator Ht-, but the transformation matrix St does 
depend on time. Opposite situation takes place in the case of the second solution X = t, 
i.e., diagonal form Hf is time-dependent and transformation matrix 5*,- do not depend on 
time. 



IV. SOLUTION OF THE RICCATI EQUATION. 



Recently, the solution of the model specified by the Eq. (1) was given in the case when 



N 



T-L = I^C^ and the operator He, V are defined as (for detailed discussion see^) 



n=l 



N 



N 



He = ^uJnCr3{n) 

n=l 



(25) 



n=l 
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where Un and gn are certain constants defining the frequencies and the coupling constant of 
the spin-bath, respectively. For every n < N the operator a^{n) is understood as o"3(?7,) = 
J2 ® ... ® o"3 (8> ... ® where J2 is a 2 x 2 identity matrix and (T3 is the usual Pauli matrix. 
Note that the operators above commute, i.e., [He., V] = 0. We will extend results obtained 
in paper ^ to the arbitrary operators that commute with each other. 

If one assumes that the operators He and V commute, then they have common set of 
eigenvectors. Henceforward, we will assume that the eigenvalues of He and V are all discrete 
and not degenerated. Furthermore, the spectrum of a given operator A will be denoted by 
(J {A). As a summary, one can write the following eigenvalue problems for He and V: 

HE\<Pn)=E^\<Pn), V\<Pn) = Vn\(t>n)- (26) 

where En G ct^He) and Vn € c^V). Here, the index n goes through the set of all integer 
numbers or through each subset of it. We wish to emphasize that the assumption of the 
discreetness and no degeneration of the spectrum of the operators in question is not crucial 
in our analysis. In fact, it can by easily overcome (however, we will not address this technical 
issue in the current manuscript). 

Note, the solution X of the Eq. (4) is a function of the operators H± defined in the 
Eq. (4). Since the operators He and V commute with each other, so are the operators 
Therefore, [X, H^] = and the Riccati equation (4) can be simplified to the following form 

aX^ + 2VpX - a = 0. (27) 

One can observe that if a = then Xq = 0^; is the solution of the Eq. (27). No wonder since 
in this case matrix Hqe is already in the diagonal form. Yet, the operator Xq may not be 
the only solution of the equation Rhqe [^] = 0. 

In order to obtain the solution for the case where a 7^ we will apply the spectral theorem 
for Hermitian operators. From the Eq. (27) one can readily see that X = /(Vg), where the 
function / takes the form 

a 

and (T{Vp) = a{V + /3Ie). One may also write the operator f{V) in the equivalent way: 
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fiv)= Yl fim(M- (29) 

Thus, in view of the Eq. (26) and (29) the X takes the final form 

X = 5^ x„ I (/)„)(</.„ I, (30) 

n 

where abbreviation x„ = f{Vn + /?) was introduced. Note that x„ are the eigenvalues of the 
operator X, that is X\(j)n) = Xn\<f)n)- The operator above is Hermitian, i.e., X = X"^. Since 
for every parameter a and (3 the function / specified by the Eq. (28) takes positive values, 
thus the eigenvalues Xn are positive. Therefore, the solution we obtained in the Eq. (30) 
is positively defined operator. Moreover, as was pointed out in^ there exists at least one 
more solution of the Eq. (27), in the case where He and V were chosen to be the ones that 
describe the spin-bath (see Eq. (25)). This situation also occurs when the operators He and 
V have more general form, then the second solution is given by 

/(A) = ^ , A G a{Vp). (31) 

a 

It can be easily verified that this function determins the negatively defined operator, that is 
to say 

X = ^x„|0„)(0„|, (32) 

n 

and Xn '■= f{Vn + /?). Observe that X = X'^ , i.e., this operator is also Hermitian. Note also 
that /(A)/(A) = —1, for A G o"(Va). This may be verified directly, or may also be thought 
of as a consequence of the Viet 'a formulas, if one thinks of / and / as the solution of the 
quadratic equation: ax^ + Ax — a = 0, for A G (^(Vs). We want to emphasize however, that 
the Riccati Eq. (27) is not a simple quadratic polynomial, even though it may appear so. 
As a result there might exist other solutions of this equation that are not specified by the 
well-known formula for roots of the quadratic equation. 

At this point natural questions may by asked, for instance which operator, X or X 
(or any other, if it exist) should by used to diagonalize the block operator matrix Hqe'^ 
What difference (if any) does it make? Of course, each solution of the Riccati equation will 
diagonalize the operator matrix with which it is associated. Nevertheless, in certain cases 
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it may by convenient to chose one solution instead of the other. For example by studying 
limiting cases like a — ?■ 0. Indeed, if /3 7^ 0, then we find that (recall that / ~ 1//) 



lim / = 0, lim / = —00. 

a— >-0 a— >-0 



(33) 



Therefore, X ^ Oe while X — )■ —00 as a goes to 0. This means that the first solution is a 
continuous function of the parameter a, including the a = value even though the second 
operator, i.e., X does not exist in that point at all. As a result, if in the process of analysis 
one decides to use the second solution X then one might meet serious problems taking the 
limits a 0. Furthermore, as was mentioned earlier the first solution is a positively defined 
operator, thus it is more suitable to manage. Henceforward, we will restrict further analysis 
to the operator X given by the Eq. (30) only. 

V. THE EXACT REDUCED DYNAMICS 

We now use the solution above to construct the exact reduced dynamics of the model 
described by Hqe. Obtaining the exact reduced dynamics of the model under consideration, 
namely the one defined by the Eq. (1), can be easily accomplished using the Eq. (3), as we 
mentioned earlier. We begin with the construction of the evolution operator Ut generated 
by the Hamiltonian Hqe- Foremost, let us recall that the evolution operator Ut may be 
computed by applying the following formula 



where Hd stands for the diagonal form of the Hamiltonian Hqe and in agreement with 
Eq. (8) it takes the form 



Ut = Sxexp{-iHdt)S 



(34) 



H+ + aX 



Oe 



(35) 



H_-aX 



where X is given by (30). Next let us observe that the inverse operator is 



S 



1-1 _ 



g{X) Oe 
Oe g{X) 



x 



(36) 



G{X)S[ 
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where g{X) is the function of X and it is given by 



9W 



(37) 



1 + A2^ 

Note that [X, V] = 0, thus X and V have the same eigenstate, namely Since X = /(Vg), 
the eigenvalues Xn of X are given by x„ = f{Vn + (3). Obviously, <j{X) = f{a(yi3)). The 
same arguments lead to the conclusion that g{X)\(f)n) = {g ° f){En + this implies 

that a{g{X)) = {g o /)(cr(Vg)). Using the Eqs. (34) - (36) one may finally write the form of 
the evolution operator of the total system, it reads 



Ut = G{X) 



Ut+X'Uf {Ut-Ui)X 

{ut-ur)x u^ + x^ut 



(38) 



where := exp(— 2(if-|- ± aX)t). Equation (38) represents the evolution operator of the 
total system Q + E and it can be easily applied to any (not only to the factorable ones) 
initial state pqe of that system, since it is written in 2 x 2 block operator matrix. At this 
point reduced dynamics PQ{t) may by obtained, it has the following form 



PQ{t) = TiEiUtpQEUj), (39) 

where Tr^; is the partial trace operation. Note that in general case this may not be easy 
to accomplish, even though the Eq. (39) might indicate so. The reason for that is that to 
compute partial trace Tr^; one needs to write operator pQE{t) in the 2x2 block operator 
form. 



VI. KRAUS REPRESENTATION 

It is interesting to see how the results of the previous section are related to the standard 
description of the completely positive map via so called Kraus sum representation: 

PQit) = Y,K,{t)pQK,{t)\ (40) 

where the Kraus matrices K^{t) satisfy following completeness relation '^^K^{t)Kfj^{ty = 
Iq. It is well established that it is possible to derive the Eq. (40) from the Eq. (39) assuming 
that no correlations between the system and its environment are present initially -'^\ The 

12 



generalization to the case when initial state is not factorable is also possible^'' (see discussion 
below). Nevertheless, in practice finding the Kraus matrices is impossible in most cases. We 
will show how to construct those matrices for the system we study. To accomplish this goal 
let us rewrite Sx in the following manner 




® \4>n){4>n\ 

(41) 

n 

We also used resolution of the identity Ie in the basis, that is to say Ie = J2n I0n)(0n|- 
Note the inverse operator S^^ can by written in a similar fashion, namely 

Sx' = J2^n'®\<l^n){<l^nl (42) 
n 

where is the inverse matrix of F„. Since we have det(F„) = 1 + > it always exists. 
Due to the fact that F'^ = F^/det(F„) holds we can rescale F„ namely F„ 1 / a/ det(F„)F„ 
so it becomes the unitary operator. Furthermore, the Hamiltonian (35) may be rewritten as 

n 

where := {E^ ± (3) ± aXn and := En ± Vn are the eigenvalues of (note cr{H^) = 
a{HE ±V))). Combining Eq. (34) and Eqs. (41) - (43) we obtain 

Ut = Y,Un{t)^\<Pn){<Pn\, (44) 

n 

with the unitary matrices Unit) = exp{—iHnt) and Hn = F~^HnFn. One can see from the 
form of the evolution operator above and Eq. (39) that we finally have 

PQ{t) = Y.PnUn{t)pQUn{t)^ 

(45) 
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In the Eq. (45) p„ := (0n|ps|0n) and the Kraus matrices are defined as Kn(t) := y/p^.Un(t). 
Therefore, the operator sum representation of the model we study is found. 

A. Connection with the Riccati diagonalization 

Interestingly, the explicit form of the matrix Hn reads 



and it does not depend on the eigenvalues x„ of the operator X . Therefore, one may draw 
the conclusion that the Kraus matrices obtained in Eq. (45) do not depend on To solve 
this puzzle, notice first that in order to compute the Krause matrices one needs to determine 
the "evolution" operator f/„(t). The later requires diagonalization of its "generator" Hn. 
Because Hn = F'^H^Fn the dependence of the Kraus matrices on the eigenvalues x„ is "hid- 
den" in a way, in the similarity matrix F„. It is important to realize that the diagonalization 
procedure Hn Hn differs from the standard diagonalization routine, which is based on 
the eigenvalue problem for the operator Hn- This new kind of algorithm is called Riccati 
diagonalization and was recently introduced in"". One may see that it arises naturally in our 
analysis. Let us also recall that the similarity matrix F„ is composed of the nth eigenvalue 
Xn of the operator X, which is the solution of the Riccati Eq. (27), while the eigenvalues of 
the matrix Hn are the solution of the following characteristic equation associated with H^- 



We have already found its solution indirectly in the Eq. (43). The roots of this equation are 
given by A.^^ = hn- The corresponding eigenvector (not normalized) reads 



Therefore, the similarity (not unitary) matrix G„ = (A^, A^) the one that also diagonalizes 
H„ takes the form 




(46) 



- ATr(i7„) + det{Hn) = 0, 



(47) 




(48) 
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G„ = I '''^ ""^^ I . (49) 

The connection between similarity matrices and is following. If F„ = (^Ji,^^), then 
^n — ^n- -^^ and = A^, since = —1. Therefore the matrix is also composed 

of the eigenvectors of yet the matrices Gn and F„ are not similar, ie., the invertible 
matrix P„ such that Gn = PnFnPn^ does not exist. It immediately follows from the fact 
that Tr(G'„) ^ Tr(F„) as well as det(G„) 7^ det(F„). Since the similarity transformation 
preserves the trace and determinant, thus the matrices G„ and F„ can not be similar. From 
the consideration above one may easily grasp the main difference between the two methods. 
In the standard approach one needs to determine the eigenvalues of if„ as was mentioned 
above {i.e., one need to solve Eq. (47)). On the other hand to use the Riccati diagonalization 
schema to our advantage we need to find x„ {i.e., resolve the Riccati Eq. (27)). Note that 
in the case of the ordinary matrix, the Eq. (27) becomes the quadratic equation, but it 
differs from the characteristic Eq. (47). In our model the relation between A^ and x„ may 
be summarized as 

A± = (E± ± /3) ± axn. (50) 

B. Initial correlation 

Let us notice that if the initial correlations between system of interest and its environment 
are present, i.e., pqe takes the form 

Pqe = ^lijPQ^ Pe, (51) 

ij 

for some not factorable complex number 7jj, namely 7jj 7^ Yili^ then reduced dynamics (39) 
cannot be written in the form (40). This not so obvious, since there are cases in which, 
even though initial correlations are present the reduced dynamics can still be written in the 
Kraus form'\ 

Yet, for a finite dimensional environment there exists a simple criterion'. It allows one 
to verify when the given state PQ{t) possesses the operator sum representation (40) if initial 
correlations are present. The necessary and sufficient condition for the later to holds true 
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for any initially correlated state is that the joint dynamics has to be locally unitary, i.e., 
the evolution operator Ut for the total system needs to be of the form 



Ut = UQ{t)(E)UE{t). 



(52) 



The operator UQ{t) {UE{t)) describes the evolution of the system Q {E) alone. Observe 
that if the evolution of the total system is not locally unitary that this not necessary need 
to implies that pqit) does not posses Kraus representation for particular initially correlated 
state pQ. From the Eq. (44) one can readily see that the evolution operator (38) does not 
have the form (52) as one may expected. Note that from (51) and Eq. (44) we obtain 



where e^- = (0^1 p-^| Therefore, even if in this general case the reduced dynamics can 
not be written in the operator sum representation, one can still describe the evolution of 
the system in a manageable way. 



In this paper we investigated the Riccati algebraic equation associated with the Hamil- 
tonian defining the process of decoherence in the case of one qubit. It was shown that if the 
environment is r— symmetric, where r is antilinear involution, then r is the solution of the 
Riccati equation under consideration. We indicated that even though the solution of the 
Riccati equation has been found it can not be applied to obtain the reduced dynamics, due 
to the problem with standard procedure allowing one to solve Schrodinger equation. We 
wish to emphasize that this result does not contradict with the previous paper \ where we 
claim that the solution of the Riccati equation enables one to rewrite the evolution operator 
generated by the Hamiltonian Hqe as 2 x 2 block operator matrix. Of course, the reason 
that problems occur is that the solution of the Riccati Eq. is antilinear. 

Furthermore, we provided a full resolution of the problems introduced in^ for the case 
when operators defining the environment commute with each other (regardless of the exis- 
tence of any symmetry in the system). We also derived the operator sum representation for 
that model assuming no correlations between the systems are present initially. Moreover, 




(53) 



VII. SUMMARY 
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we also showed how to obtain the solution if the initial state of the total system in not fac- 
torable. This result is a direct generalization of the system discussed in^. We also indicated 
that the recently derived schema of so called Riccati diagonalization arises naturally in the 
model we considered. 
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